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In this paper, we present a method for the numerical minimization 
of the Mumford-Shah functional that is based on the idea of topological 
asymptotic expansions. The basic idea is to cover the expected edge 
set with balls of radius e > and use the number of balls, multiplied 
with 2e, as an estimate for the length of the edge set. We introduce a 
functional based on this idea and prove that it converges in the sense 
of F-limits to the Mumford-Shah functional. Moreover, we show that 
ideas from topological asymptotic analysis can be used for determining 
where to position the balls covering the edge set. The results of the 
proposed method are presented by means of two numerical examples and 
compared with the results of the classical approximation due to Ambrosio 
and Tortorelli. 

1 Introduction 

Let ft he a Lipschitz bounded open domain in M^. We assume that a possibly 
noisy image on fl is given, represented by a real-valued, bounded function / on 
ri, whose values /(x), x £ Q, correspond to the intensity of / at the pixel x. 
In order to segment, and denoise at the same time, the image /, Mumford and 
Shah jSBj introduced a variational model, which is based on the assumptions 
that the different "objects" in the image give rise to homogeneous regions that 
are separated by the objects' projected silhouettes. Moreover, these silhouettes 
in general correspond to discontinuities in the image /. By this reasoning, they 
proposed to minimize the functional 
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[ \S/u\^dx + n\K), (1.1) 



taking as variables the function u E L'^{n), the denoised image, and the com- 
pact set if C ri, the set of edges or silhouettes. Here, 'H}{K) denotes the 
one-dimensional Hausdorff measure of the set K] in the case if is a regular 
(rectifiable) one-dimensional set, this is precisely its length. The parameters a 
and /3 that appear in (1.1 ) are positive constants determining the weight that is 
put on the regularity of the denoised image u and the length of the edge set K. 

In order to show the existence of minimizers of the Mumford-Shah func- 
tional, a weak formulation depending only on one variable has been introduced 
by De Giorgi et al. (26^. In their model, the set K is replaced by the jump 
set Su of the function u E SBV{il), the space of special function of bounded 
variation on D,. Still, this reformulation of the functional provides no method 
for the actual numerical computation of minimizers. Thus, various approxi- 
mations of the functional F have been proposed, most of them based on the 
concept of F-convergence. Ambrosio and Tortorelli 3 proposed a variational 
model in which they replace the set by a continuous function v with values 
close to near Su, and values close to 1 away from Su- For other approxima- 
tions of the Mumford-Shah functional in the sense of F-convergence, we refer 
to Braides et al. 13 , where the authors propose approximations by a family of 
non-local functionals. Approximations by finite-difference schemes, inspired by 
the original discrete model of Blake and Zisserman [TD] , have been considered 
by ChamboUe (ini [U] , and by finite-elements schemes by ChamboUe and Dal 
Maso [IH] or recently by Aubert et al. [S] . In the work by Koepfler et al. [35] and 
Dal Maso et al. I25J , region growing and merging methods have been proposed. 
For detailed analysis of the Mumford-Shah model, we refer reader to the book 
by Morel and Solimini [35] . We also refer to the books [3 [39l HQ] , where some 
of the above mentioned results are shortly discussed. 

In the following, we will show how the problem of minimizing the Mumford- 
Shah functional F can be approached by means of an asymptotic expansion 
that is based on the idea of topological asymptotic analysis. In its original for- 
mulation (see, e.g., Sokolowski et al. [H], Garreau et al. [SD], Feijoo et al. [^). 
this theory investigates a variation of some objective functional depending on 
subsets of K" with respect to the subtraction of a small ball from the consid- 
ered domain. This variation is a scalar function, called the topological gradient 
or the topological derivative, and its largest negative values indicate positions, 
where it is good to remove a small ball. In [S] , Amstutz proposed to modify the 
definition of topological gradient and provide variations of a given functional 
with respect to change of certain material properties, and not a domain topol- 
ogy. Recently, topological asymptotic analysis has been also applied for image 
processing purposes by Auroux et al. [H] |9] and by Muszkieta 07] . 

In order to apply the theory of topological asymptotic expansions, it is nec- 
essary that the functional to be minimized depends solely on the set K. This 
can be achieved in the case of the Mumford-Shah functional by noting that 
the minimizing pair (uq,Ko) is uniquely determined by either of the compo- 
nents uq and Kq: The set Kq coincides with the jump set Suq of the function 
uq. Conversely, uq can be computed from Kq by solving the partial differential 
equation 

u — aAu — f in , 

du (1-2) 
^^=0 md{n\Ko). 
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Now consider the functional 



JiK)^F{u{K),K), 



where u{K) denotes the solution of (1.2) with Kq replaced by K. Then the pair 



{u{Ko),Kq) minimizes the functional F, if and only if J{Kq) < J{K) for all 

K c n. 

The idea is now to use a gradient descent like approach to the minimization 
of the functional J. Starting from an initial guess K of the edge set (for instance 
K = one adds to K those points, whose inclusion would lead to a near to 
maximal decrease of the cost functional J. More precisely, one adjoins to the 
set K small balls of radius e > centered at the points x G \ if and tries to 
compute an asymptotic expansion of the form 

J{K U B^{x)) - J{K) = c{e)GK{x) + o(c(e)) 

for some functions c: M>o — > M>o and Gk : \ if — > M. Those x G \ if for 
which Gk attains the largest negative values are then added to the set K. This 
process is iterated, until the function Gk becomes non-negative everywhere. 

In the case of the Mumford-Shah functional, this procedure cannot applied 
directly, because the functional is infinite whenever K contains a set of positive 
Lebesgue measure. We therefore propose to use a different, though related, 
functional for the computation of the asymptotic expansion, which is based 
on an approximation of the one-dimensional Hausdorff measure: The number 
of balls of radius £ > that are required to cover a set if, multiplied with 2e, 
provides a good estimate of H^{K) for e small enough and if sufficiently regular. 
In the following, we introduce this approximating functional Jg.^. However, 
because we later prove the F-convergence of this functional to F, it is necessary 
to let Jg i^ depend on two functions, the function u and a piecewise constant 
edge indicator function v: 

For each finite set y C and each < k < 1 we define the function 
vy,k : -> M by 

VY,Kix) ■■= { 



if a;eUaey^e(2/), 
else. 



For every v G L'^{n) we define 

m,,,(«) := M{n°(Y) ■.YcR', vy,^} . 

Here we set m^^i^{v) := +oo, if u 7^ vy.k. for any Y C M^. Note that, for a given 
function v, there might exist different sets Y CM.'^ such that v = vy.K- 

Finally, we introduce the family of functionals Je,K : L'^{n) x L'^{n) — > M, 
defined by 

J,,,{u,v) -.= 1 [ {u-f)^dx+^ [ v\Vu\^dx + 2(5em,,^{v) (1.3) 

if M G iJ^(ri), and Jg{u,v) :— +00 otherwise. 

In Section [2[ we derive an approximation of the functional J^^^ that allows 
us to compute an approximation of the minimizer using the idea of topological 
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asymptotic expansions. We show that 

1 — K, 

Je,K{u,VY\j{y},f,) - uy,^) « -e^Tra \Vu{y)\^ + 2e/3 . 

(cf. Theorem [2l|). 

In Section [3j we show that the functional J^^k is indeed an approximation of 
the Mumford-Shah functional in the sense of F-convergence. More precisely, if 
we choose k = K,{e) in such a way that K,{e) = o(e) as e — >■ 0, then 

F = T- lini J^.^(^) . 

In particular, this implies that the minimizers of Je,K{e) converge to a minimizer 
of as £ tends to zero. The adopted proof is based on the methods used 
for proving the F-convergence of the Ambrosio-Tortorelli approximation of F 
(see OE]). Finally, in the Section [ij we propose an algorithm to minimize 
the functional Je,K- We compare numerical results obtained with this algorithm 
with results obtained by minimization of the Ambrosio and Tortorelli model [3] . 



2 Topological Asymptotic Analysis 

In this section, we derive the topological asymptotic expansion of the functional 



Je,K, defined in (1.3). This expansion will allow us to find its approximate min- 
imizers. We assume that J7 C is a bounded and open Lipschitz domain and 
that / e H^(n) is given. We define the functional G : L^{n) x L^{n) ^ [0, +oo], 

^ Jn ^ Jn 

Now assume that K is an open subset of f2 and < /t < 1 satisfying an < 1. 
We define the function v: il M by 



v{x) 



K ii X £ K, 
1 else. 



Using standard methods of variational calculus one can show that the map- 
ping u I— > G{u,v) attains a unique minimizer in H^(fl), which we denote by 
u. 

The main result of this section is stated in the following theorem: 



Theorem 2.1 Let K C fl. For y E n \ K and e > define the functions 

v^y^'\x) 



K ifxeKUB^iy), 
1 else, 



and 

._ argminG(fi,v(*'^)). 

ueH^{n) 

Then for all compact subsets L d Q \ K we have 



lim sup — 



1 -I- K 



= 
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For the remaining part of this section we assume that the compact set L C 
\ ^ is fixed. Moreover we define 

6 dist(L, dnuK)=^ min{ \x-y\:xeL, y ednuR] . 
3 3 

In addition, we consider the set 

L:^L + BsiO) ■.= {x + y:xeL, \y\ < S} . 



Before we give the proof of the Theorem |2.1[ we need to introduce some 
auxiliary results. First, we recall that the assumptions that u and u^y'^^^ are 
minimizers of G{-,v) and G(-, u'^''^-'), respectively, imply that 



{u — f)ip + aviyu, V(/3) dx — , 



(2.1) 



for all ip e H\n). 

We first we need a regularity result for the function u. 

Lemma 2.1 The function u satisfies 

ueCUn\K). 

Moreover, there exists a constant Ci only depending on L and K such that 

^ C'i||/||ffi(o) . 

Proof. Because by assumption / e H^{il) and the function v is constant 
on il\ K, it follows from standard theorems on the regularity of solutions of 
elliptic equations that u G _ffj^^(ri \ K) (see [3TJ Thm. 8.10]). Moreover, there 
exists a constant C only depending on L (and therefore on L and K) such that 
||w||^3(£) < C'||/I1hi(o)- Thus the Sobolev embedding theorem [TJ Thm. 5.4] 
implies the assertion. □ 

As a second step, we need iJ^-norm and L^-norm estimates of the difference 
_ First we show that the 77^-norm of the difference u^^'^^ — u on the 
whole domain is of order e. 

Lemma 2.2 There exists a constant C2 > only depending on L, k, fl, and 
K such that for all e > with e < S and y (z L the estimate 

\\u^^''^ -u\\HHn)<C2\\f\\HHn)e 

holds. 



Proof. Computing the difference between the two equations in (2.1) and 
using the definition of v'-y'^\ we obtain that 

(it(*'^) -ii)(^da; + a / z;'*'''' (V(u(*'''' - u), V(^) 



a{l - k) (Vu, V</?) dx 
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for all ip G }i^(SX). In particular, it follows with Lp — u^^^'^^ — u that 



aK\\u^y^^^ - M||^i(m < a(l - k) \ {Vu, V(u(^^^) - u)) dx 

< a{l - K)||Vu||L2(s^(f;))||w(^''') - u|lHi(a) • 
Moreover Lemma [O] implies that 

l|Vu||i2(B^(jj)) < ||Vu||ioe(B^(g))V^e < \\Vu\\^^f^i^y^e < Ciy/^\\f\\H^n)e . 

(2.2) 

Setting C2 := (1 — k)\/ttCi/ k, the assertion follows. □ 

Lemma 2.3 There exists a constant C3 only depending on L, k, fl, K , and 
||/||ffi(n) such that 

for every y € L and < e < 5. 

Proof. Let g e H^{fl) satisfy ||<?||hi(si)_^ 1 and assume that lo^^''^) and w are 



the unique solutions to the equations (2.1) but with given g instead of /, that 
is 

(2.3) 

{w — g)(p + aviyw, "S/ip) dx — Q , 

JQ. 

for ah <y9 e H^{VL). 



Taking ip — u;'^^''^) and ip = w in the first and the second equation in (2.1 1, 



respecti vely, and next subtracting these equations from the corresponding equa- 



tions in (2.3 1 with ip = u^'^''^^ and if = u, we obtain 



w^y-'^\fdx- / u^^''') 5 = , 



w f dx — I ugdx = . 
n Jn 

In particular, 

/ {u^y-^^ -u)gdx^ [ {w^y-"'^ - w) f dx . 

Next, we note that 



(u(*'^)-u)gda;= i ( / {u^y^^^ ^u)gdx+ I (w'-y^"' - w) f dx 
n Jn 



2 

1 I (^(u«'^)+u;«'^))-(w + u;)) {f + g)dx 



^ I {u^y^^^ ~u)fdx-- I {w^y^^'> ~w)gdx. (2.4) 
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if — u^'^''^\ respectively, we obtain 



Computing the difference of the two equations in (2.1 1 with ip = u and 



a(l — k) 



AppHcation of the Cauchy-Schwarz inequahty to the above formula, and next, 
the estimate (2.2) and Lemma 2.2 yields the inequality 



< a(l - k) [\\yu\\L2(B,iv))\\ 

< a(l - {VnCiC2 + ttC!) \\f\\j,.^n^e 



""11^1(^2) 



ll^^llL^(i3.(y)) 



(2.5) 



In a similar manner, using the assumption that ||5|lHi(a) = 1, we can show that 



(w 



w 



)gdx < a(l - k) (V7FC1C2 + nCf) , 



(2.6) 



and 



(u + w)) (/ + g) dx\ 
< ail ~ k) (V^CiC2 + i^Cl) 2(||/||?,,(,,) + l)e^ . (2.7) 



Finally, combining ([23]), ( pi3| and ( p?7| ) with ( [2^ , we obtain 



Therefore, we have 

u\\H~^n) ■= sup< 



llJS.e) 



{u 



u)gdx -.geH^in), \\g\\ 



< 1 



Lffi(f2) 

<c(||/||?^i(^,) + l)e' (2.8) 



with c = 3a(l — k) {\pi^C\Ci + ttC^) /2. Now, estimates from the theory of 
Hilbert scales (see [33l Thm. 9.4]), Lemma 2.2 and (2.8) imply that 

<c{\\f\\lHn) + ^yCMm(n)e. 



Therefore, the desired estimate holds with C3 — ^/cdl/H^i^jj^ + l)C2\\f\\H^{n)- 
□ 
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In the next step, we need to derive an estimate for the function u^^'^) on 
the boundary dBs{y). To do this, we follow Vogelius et al. 03], where such 
an estimate has been derived for the solution of the homogeneous Helmholtz 
equation with Dirichlct boundary conditions. 

We first introduce the Green function corresponding to the equation u — 
aAu = f on ^ \ K with Neumann boundary conditions. That is, the function 
N{-,y) solves the problem 

' N{x,y)-aA^N{x,y)=Sy{x) x&n\K, 

^^{x,y)^0 x€din\K), 

for y E H.. We note that N{-,y) can be written as the sum of the fundamental 
solution corresponding to the equation u — aAu = dy, denoted by T{-,y), and 
a corrector function h{-,y), which is chosen in such a way that the normal 
derivative of N{-,y) vanishes on the boundary d{il \ K). The function r(-, y) is 
given by 

for all X, y gM.'^, such that x ^ y. Here Kq denotes the modified Bessel function 
of the second kind (see, e.g., [BU p. 490]). Furthermore, the function /-Co has an 
asymptotic expansion of the form 

Ko{z) = -lnz + ln2-7 + 0(z2|lnz|) 

for z — 0, where 7 denotes the Euler-Mascheroni constant (see, e.g., [121 
Ch. 51]). Therefore, we conclude that r(-,?/) can be approximated as 

T{x,y)^^(-\n{\x-y\) + Una + \n2-^)+0[\x-y\^Mx-y\\), (2.9) 

when I2: — j/| — > 0. Moreover, we observe that T{-,y) has the same singular 
behavior as the fundamental solution of the Laplace equation 

*(a;,2/) = M\x - y\) 



defined for all x, y e M , such that x ^ y. We need the approximation (2.9) 
in order to be able to apply standard methods of potential theory to derive an 
estimation for the function m^^'^) on dB^{ij). Such way of proceeding is common 
when dealing with problems of this kind (see, e.g., Colton and Kress pPl I21j ). 

Lemma 2.4 There exists a constant C4 only depending on L, K , and VL, such 
that for every point y E L. < e < S, and y G f2 satisfying e < \y — y\ < 2e the 
estimate 

\\Ni;y)\\mB^)<Cis\lne\ 

holds. 

Proof. Since by elliptic regularity h e C°° {Q \ K,i} \ K) and 

dist{y,n\K) > 3(5 > 3e, 
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we have 

\\hi-,y)\\L^B,{y)) < v^\\H-,y)\\L-°(n\K)£ < y^cie 

with ci > only depending on L, K, and fl. The estimate \\^{-,y)\\L2{Bc{y)) ^ 
C2e|lne| for some C2 only depending on e and for all y satisfying e < |y — y| < 
2e can be easily derived in the polar coordinate system. Therefore, from the 



Minkowski inequality and (2.9) we get that 



\\N{-,y)\\L^B,(y)) < \\h{-,y)\\L^B,{y)) + \\T{-,y)\\mB,iy)) 

< y/ncie + C2e|ln£| < C4e|lne| 
with C4 chosen slightly larger than C2- □ 

Lemma 2.5 There exists a constant C5 only depending on L, K , VL k, a, and j , 
such that for every point y G L, < e < S, and y G ^ satisfying e < \y — y\ < 2e 
the estimate 

u^y^'^y) - uiy) - ail ~ K.) f u^y''\x)^{x,y) ds{x) <C^e^/'^ 

JdB,{y) on 

holds. 

Proof. Computing the difference of the integral representation formulas of the 
functions u and u(*'^\ we obtain 



u^y''\y)-u{y) 



9B,{y) 



r~ . dN du^y^^'^ \ 



(2.10) 

a / N{x,y)^{u'^y'^^^{x)~u+{x))ds{x) 



dK 



dn 

-a I N{x,y) f{x)dx . 

for all y e n\{KU B^iy)). 

Using the transmission condition that u^^*'^) satisfies on dB^{y), the Green 
formula, and that 

N{x,y)~al^N{x,y)^{) (2.11) 
for X £ Bg{y) and y G Vl\[K U B^{y)), we have 



a I N{x,y) — {x)ds{x)^aK N{x,y) — (x) ds{x) 

dB,(y) dn jQB,(y) dn 



= an N{x,y)Au^y^'''>{x)dx + aK N{x,y),Vu^y'''> {x)) dx 

JB^iy) JB,(y) 

= {1-k) [ N{x,y)u^y^'^x)dx~ [ N{x,y)f{x)dx 

JB,(y) JBeiy) 

+ aK [ u''y''^Ux)^^{x,y) ds(x) . 
JdB^iy) dn 



(2.12) 
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Using similar arguments as above and applying the Cauchy-Schwarz inequality, 
we obtain 

-a f N{x,y)-^{u'-y'''^^{x) -u+{x))ds{x) 
JdK 9n 

= —an / N (x , y) — iu"^ '^"^ (x) — (x)) ds(x) 
JdK on 

= ail-K,) [ {u^y^'\x)-u{x))N{x,y)dx ^^-^^^ 
Jk 

< a(l - K)\\u^y''^ - u\\L2^K)\\N{-,y)\\L2,^K) ■ 
Because dist(y,_ft') > dist(y,iir) — 2e > S, there exists some c > such that 



W^i' J y)\\L^{K) < c. Application of Lemma 2.3 therefore yields 

|!««'^) - u\\LHK)\\N{-,y)\\mK) < c\\u^^''^ - u\\L2^n) < C^ce-I^ . (2.14) 



Taking into account ( |2.12[ ), ( |2.13| ), and (|2.14| in (|2.10|), we obtain 



u^^''\x)^{x,y) ds{x) 



dBAy) 



dn 



-{1~k) / u(y''\x)N{x,y)dx 

+ {l~a) I f{x)N{x, y) dx + 0{e^''') , 
•IB-Ay) 

the last term being bounded by a{l — k)C^c£^/'^ . To estimate the remaining 
integrals, we note that 

Ik^^'^^IU^Csas)) < v^ll/IU-(0)£. (2.15) 
Next, using the Cauchy-Schwarz inequality and Lemma|2.4| we get 



u(^'^)(x)7V(x,y)dx< h»'^)|U2(B.(.5))||A^(-,y)|U2(B,te)) 

BAv) 

< V7r||/|U-(o)C4e^|lne|. 

In a similar way, we show that 

f{x)N{x,y)dx < \\f\\L^BAv))\\^i-^y)\\LHBAyy) 

< V7r||/||L-(o)C4e^|lne| . 



Be{y) 



Therefore, we obtain the desired estimate with C5 slightly larger than a{l — 

k)C3C. □ 



Lemma 2.6 There exists a constant Cq only depending on L, K , k, a, and j , 
such that for every point y €z L, < e < S, and y ^ fl satisfying e < \y — y\ < 2e 
the estimate 



uiV'^) (x) ^ {x, y) ds{x) - I u^y'-^ (x) ^ {x, y) ds{x) 
dBA-y) "'^ JdBAy) "'^ 



< C6£^|lne| 



holds. 
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Proof. Denote H{x,y) := h {x,y) +V{x,y) - <^{x,y). Then N{x,y) = 
^{x,y) + H{x,y) and, using ( |2.11[ ), we obtain 



N{x,y)^ ^H{x,y). 



(2.16) 



Application of the Green formula, (2.16), and the Cauchy-Schwarz inequality 
yields 



u^^'^\x)^^{x,y) ds{x) 



dH 



dn 



B^y) 



'^y''\x)AH{x,y)dx+ / {\/u^y^''\x),\/ H{x,y)) dx 



B^y) 



'^y''\x)N{x,y)dx+ / {Vu^y-'^\x),VH{x,y))dx 



B.{y) 



< ||w^^'^^||L2(B^(y))||iV(-,2/)||L2(i3,(y)) 

+ llVu||L2(B,(y)) ||ViJ(-,y)||i2(B^(y)) 



Using (2.15) and Lemmas 2.2 and 2.4 it follows that 



dB,(y) 



Sy-'\x)"^{x,y)ds{x) 



< 



C^4V^||/||Lo^(i2)e'|ln£| + (C2+CiV^)|l/||«i(0)||Vff(-,y)|U2(B,(g))e (2.17) 



Using the explicit form of r(-, y) and <i>(a;, y), and next the Taylor expansion of 
the function Kq (see T2, Ch. 51]), we can show that 

|Vr(x,2/)-V$(a;,y)| < c£|lne| 

for some c > only depending on 5 and all x e Bg^y) and y satisfying e < 
\y~y\< 2e. Thus 

\\'^H{-,y)\\L2(B,iy)) < {ce\\ne\ + \\h\\Lo.^B,iy)))\^£ 



Thus we obtain the required estimate from (2.17) by choosing Cq slightly larger 
than C4V7r||/||L~(o)- 



□ 



From Lemma |2.5| and Lemma [2.6| and using the jump formula for the double 
layer potential (see e.g. Kress |MJ p. 67]), we have 



—a{l — k) 



■Ay) 



u^y^'\x)^{x,y) ds{x) 



< Cre' 



3/2 



(2.18) 



for y e dB,{y), with Cy > C5. 

Now we introduce the auxiliary vector valued function : M x M that 
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solves the problem 
' A0(x) = 



X e \ B(0) or a; e B(0) , 
X e dB{0) , 



(x) — K— — (x) = —Kn{x) X G dB{0) . 



(2.19) 



dn dn 
lim 0(a;) = . 

^ |a;|— foo 



Here B{0) denotes the unit ball in of a center in and n is the unit normal 
vector exterior to the boundary dB{0). The existence and uniqueness of a 



solution to the problem (2.19) can be proved using single layer potentials with 



suitably chosen densities. For details, see Ammari and Kang [3] or Cedio-Fengya 
et al. [H. 



Solving the problem (2.191 using standard methods of potential theory, we 

(2.20) 



obtain the explicit form of 0, which reads 



0(x) 



K+1 



X and <^{x) 



K + 1 \x\ 



for ah X e i3(0) and x e \ 5(0), respectively. 

The result on asymptotic expansion of the function on the boundary 

dB^{y) is stated in the following Lemma: 

Lemma 2.7 For every point y d L, < e < S, and y e dB^{y) the estimate 



u^y'^\y) - uiy) - £ ( ^ - 1 ) (0(y/e), Vii(y)) 



holds, where the constant Cy is as in (2.181 



Proof. The proof of this lemma can be proved starting from the formula ( 2.18 ) 
in the same way as in Vogelius and Volkov [43l Prop. 3]. □ 



Using all the above results, we can now prove Theorem 2.1 



2.1 Proof of Theorem \2A\ 

Using (2.1 1 with ip = u'''^^) and ip = u, we obtain that 



G{u^y'^\v'^y^''>)-G{u,v) = -^ [ f{u^y-^^ -u)dx 

2 ./o 



(2.21) 



Again using (2.1 1, it follows that 

^ ^ /(u^*-'^) - u) dx 



dx 



= -\a{l-K)l {Vu,Vu^y-'''^)dx . 

2 JB,(y) 



(2.22) 
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Next, we apply Green's formula to u G H^{^) and use the fact that u solves 
the equation u ~ f — aAu on B^{y) to obtain that 



a I {\Iu,\Ju^y^''^)dx 



'BAy) JbAv) 



{u 



BAv) 
a 



u){f-u)dx + a I (it(^'") -u)^ds 



du 



9BAv) 



dr 



BAv) 

To estimate the first integral on the right hand side of above equation, we use 



the Cauchy-Schwarz inequality and Lemma 2.3 and obtain 

{u^y'^^ -u){u- f)dx < ||m - /|U2(s^(g))||M(^'^) -u\\L2(^BAy)) 
< 2||/|Uo.(o)V^C3e^/2 . 



BAy) 



(2.23) 



Lemma 2.1 implies that Vm is Lipschitz continuous on every compact subset 
oi n \ K. Thus there exists some constant c such that 



sup |V^(^)-V"(^)I<, 
x€BAy) F-y| 



(2.24) 



for all y G i and < £ < dist(L, ft \ K)/2. In particular, we have the estimate 

/ \\7u\^ dx - e^n\\/u{y)f' 
JbAv) 



{Vu - ^u{y), Vm + Vu(y)) dx 



BAy) 

< cs \\/u + S/u{y)\dx 

■J BAy) 

.3 



(2.25) 



<2cTTCi\\f\\mme' . 
The change of variable x = y+ex, application of Lemma|2.7|and the property 



(2.24) yields 



du 



9BAv) 



(u^y-^>(x)-u(x))—(x) ds(x) 
' on 



du 



e I (u'^y^'^^ (y + ex) — u(y + ex)) t— (y + ex) ds(x) 
ldB{o) on 

(- - l) Vu{yf [ m^ieS:) ds(i) + 0(e^/2) 
V'* / JdB(o) on 



(2.26) 



1 



e- I - - 1 \/u{yf / (j){i)n{£fds{x) \/u{y) + 0{e^/^) 
« / \JaB{o) J 
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where T denotes the vector transpose. Combinmg the estimates (2.23), (2.25) 
and \2.2Q\ with ( |2.21[ ) and ( |2.22[ ), we obtain 

.1 ,^ . fl 



-£2^a(l-«)7r|V7.(y)p + 0(£5/2), 



dB 



x)n{i) ds{x) j Vu(y) (2.27) 



Inserting the exphcit formula for the function (/), given in (2.20), to the 

1 - K 



asymptotic expansion (2.27), we get 



1 + K 

which ends the proof of Theorem |2.1[ 



3 The F-convergence of Je ^(e) to F 

In this section we prove that the sequence of functionals J^^k converges to the 
Mumford-Shah functional in the sense of F-convergence, if the parameter k, 
depending on e, tends sufficiently fast to zero as £ — ?> 0. 
Let us first recall the definition of the F-limit: 

Definition 3.1 Let X be a topological space and Jj : X — > [0, +oo] a sequence 
of functionals on X . Denote moreover, for x € X , by Af{x) the set of all open 
neighborhoods of x. Then the F-lower limit and the F-upper limit of Jj are the 
functionals defined as 

(F-liminf Jj )(m) :— sup liminf inf Jj{v) , 
(F-limsup Jj)(u) :— sup limsup inf Jj{v) . 

If the T-upper and lower limits coincide, we define the T-limit by 

(F-lim J,)(u) — (F-limsup Jj)(m) = (F-liminf Jj)(m) . 
J" j j 

In metric spaces, the F-limit of a sequence of functionals can be characterized 
by means of the following result: 

Lemma 3.1 Let X be a metric space and Jj' X — !■ [0, +oo] a sequence of 
functionals on X. Let moreover J: X — )■ [0, -|-cx)] and let X be a dense subset 
of {u € X : J{u) < oo}. Assume moreover that for every u £ X there exists a 
sequence Uj — )■ u with Uj G X such that J{uj) J{u). Then J = F-limj Jj, if 
the following conditions hold: 

1. For every u £ X and every sequence Uj — ^ u with Uj £ X we have 

J(u) < lim inf Jj (mj) . (3.1) 
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2. For every u G X and every S > there exists a sequence Uj — > u with 
Uj G X for all j such that 



J{u) > limswp Jj{uj) — 5 



(3.2) 



Proof. The result follows by combining the standard metric characterization 
of the F-hmit (see, e.g., [211 Prop. 8.1]) with a diagonal sequence argument. □ 



Define 



F{u,v) 



(u - ff dx + 



\Vu\^ dx + pn^{Su) 



(3.3) 



The main result of this section is the following theorem: 



Theorem 3.1 Let F and J^.n he as in (3.3) and (1.3), respectively. Assume 
moreover that K{e) — o{e) as e — >■ 0. Then we have for every sequence e.j 
that 



F = r-lim J, 



ej,K(ej) 



We now prove Theorem |3.1| using the methods introduced by Ambrosio and 
Tortorelli [2], as presented in the notes by ChamboUe [15] and the books by 
Braides [TTJ [T^]. The proof is split into three parts. In the first part, we 
will prove the lower bound, inequality (3.1), in the one-dimensional case. In 
the second part, we will extend this result to dimension 2 using the slicing 
method (see, e.g., [TTJ HI]). In the last part, we will prove the upper bound, 
inequality (3.2 ). 



3.1 The Lower Bound for n = 1 

Let il C M be open and bounded, and / € L°°(ri). We define the one- 
dimensional Mumford-Shah functional F : LF'{VI) x £^(51) -> [0, H-oo] as 

[u-ffdx+^l {u'f dx + m'^[Su) ifw=l, 

else. 




F{u,v)= {'^ J^" 2 



Because of technicalities of the proof of the two-dimensional case that result 
from the restriction of the approximating functionals Je_„ to lines, it is necessary 
to use a slightly different definition in the one-dimensional case; instead of only 
covering the edge set with balls of radius e, we also allow covers with smaller 
balls. For each finite set Y = {yi : I < i < m} of points in M we denote by 
M^^f^{Y,i^) the set of all functions w: 51 -> M for which there exist a sequence 
of positive numbers smaller than, or equal to, e, such that for all x G Q 

we have 

1 else. 



vlx) 



Furthermore we denote by 

meAv,n) ■= inf{-H°(r) : r C M, v£ M,,j,{Y,n)} . (3.4) 
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As in the two-dimensional case, we define 7fie^K{v, SI) := +00 if w ^ M^^^iY^ f2) 
for any Y cM?. 

Finally, we define the functional J^^k. ■ L'^{^) x i^(r2) — > [0, +00] as 

Je.H.{u, w) := ^ I {u- ff' dx+^ I v{uf dx + ^) ■ 

^ Jn ^ Jn 

Let {u,v) e L^(Sl) X _L^(S1), and assume that Sj — )■ 0, Uj — )■ w, and Vj — z;. 
We show in the following that 

F{u,v) < liminf . 

We may assume without loss of generality that the sequence .^^(e ) (uj, Wj) 
converges to some finite number c < +00; if c = +00, the claim trivially holds. 
Moreover we note that the continuity of the norm implies that /q(uj — f)^ dx — > 
J^{u-f)^dx. 

We first define 

T,:= {xefl:u^ H^{x- 6,x + 6) for all 5 > O} . 

In the following, we show that S is a finite set. Let therefore si < S2 < ■ ■ ■ < 
Sm G S, and let (5 > be such that the sets {si — 2S, Si + 2S) are disjoint. Because 
u ^ H^{si ~ S,Si + S) and uj — > u in L^(f2), the lower semi-continuity of the 
functional w 1— > J{w')'^ dx implies that 

liminf / (u'A'^ dx > / {u'^ dx ^ +00 . 

Js,-S Js,-S 

Thus the inequality 

limsup / {u'aY' dx essinf Vj 

2 ~ 2c 
< limsup / Vjiu'A"^ dx < — lim Jg.iuj^Vj) — — 



implies that 

essinf — >■ for all 1 < i < m . (3.5) 
Moreover, we have that 

1 ~ c 
Iimsupm5^,^(e^.)(wj,r2) < - lim Je,Mej){uj,Vj) = - . (3.6) 

Thus, for i sufficiently large, there exist sets Yj with THP{Yj) <(/?-!- l)/c such 
that Vj G Mg^. K(ej)(yj , f^)- In particular, this implies that Vj — 1 with respect to 
the L^-norm, showing that v = 1. Now the definition of the sets Mg . _K(g .)(K,-, fi) 

and (3.5) imply that for every 1 < i < m there exists some x\ G Xj with 

B^- (xp^ n {si — 6,Si + S) 7^ 0. Because Sj < 6/2 for j large enough and the sets 

{si — 25, Si + 26) are disjoint, it follows that x^'' 7^ a;|/^ for j large enough and 
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i ^ i' . Consequently, ?^''{E) < rhi,.^f^(^i,.^{vj,n) for j large enough. Because by 
construction u € H^^^{n \ E), it follows that Su C E, and thus 

H"(5„) < H"(E) < liminfm,^.„(..)(v„l]) < +c» . (3.7) 

Next we have to estimate the term Jq^j^ (^0^ dx. By assumption, the se- 
quence (wj)jgN converges to v with respect to the L^-norm. Thus Egoroff's 
Theorem (see, for instance, [281 2.3.7]) implies that for every 7 > there exists 
a set ftj C ft\Tt with C^{n \ fl^) < 7 and a subsequence Vj^. such that vj^^ 
converges uniformly to v = 1. Because Vj^, takes only the values 1 and k^-^ and 
Kg^.^ — > 0, this shows that, in fact, Vj^, (x) — 1 for x £ fly and k sufficiently large. 
Therefore the convergence of Uj^. to u implies that 




Taking the supremum with respect to 7, it follows that u e H^{fl \ Su) and 




which concludes the proof. 



3.2 The Lower Bound for n = 2 



The second part of the proof of Theorem 3.1 is concerned with showing (3.1) 
for C . The proof applies the slicing method following Braides [TTl [T2] . To 
that end it is necessary to introduce some notational preliminaries: 

We denote for every direction ^ g := {x G : \x\ = 1} by 11^ := {y G 
• (2/jO = 0} the hyperplane passing through that is orthogonal to ^. If 
A C is open, we denote by A^^^y := {t E R : y + G A} C M. the one- 
dimensional slice of A indexed by y € H^. Finally, for all w defined on 17, we 
define the one-dimensional function w^'y{t) — w{y + t£^) as the restriction of w 

to 1^5, J,. 

Next we define for every open set A C f2 a localized functional J^,k{u, v, A). 
To that end, we first localize the functional m^.^. We define 

m,^^iv,A) ■.= mi{H%Y):Y CR", vU ^ vy,.\a} ■ 
Then we define 

J,,^{u,v,A) := ^ / ff + av\\/u\^dx + 2ePTne,^iv,A) . 

^ J A 

Moreover, we define for each ^ G §^ the directional functional 

Jl^iu,v,A) := ^ J^{u- + av{^,Vu)^ dx + 2ePm,,4v,A) . 
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Finally, we consider for each ^ € y € 11^, and / C fl^^y open the one- 
dimensional functionals 

^ JI\Sa 

and 



2 _ 

where fhe^K. is as defined in ( |3.4[ ). 

We claim that for every w, -y e L'^{^)j e>0, 0<«;<1, AcO open, and 
^ S §^ the inequalities 

JeAu^v,A)>Jl,{u,v,A)>[ Jl;y{u^'y,v^'y,A^,y)dn\y) (3.8) 

hold. Indeed, the first inequality is a direct consequence of the definition of 
the involved functionals. For the second inequality, note first that, by Fubini's 
theorem, 

{u-ff + av{^,Vu)'^dx^[ ( {u^-y ~ f^'yf + av^^y{u^-y)''^ dtdH^{y) . 

Thus it remains to show that 

2£m,,,(«,A)> / 7h,,^{v^^y , A^^y) dn\y) (3.9) 

whenever v G L'^(A). In case rrii^^i^^v, A) = +oo, this inequality trivially holds. 
Else, there exists a set y = {yi, . . . , y^} C with m — m^^i^{v, A) such that 

v{x) 



He if a; e Ui=l -SeCyO : 

1 else. 



Then 

m 

2e m - ^ {{y E : B,{y,) D {y + MO ^ 0}) 

m 

> 5^ ({2/ e : B,(2/,;) n (y + MC) n A ^ 0}) 

: B,{y,) n{y + m)nA^ 0}) dn\y) . 

Moreover the definition of rh^^i^{v^'y , A^^^y) implies that 

: S,(yO n A^^y ^ 0}) > m,,,(^;^■^ A^,,) 



for all y and ^. This shows (3.9), which in turn implies (3.8). 

Now let {u,v) e L^(ri) X L^(f7), and assume that Sj — > 0, uj — ^ u, and 
Vj — 7- f . As in the one-dimensional case we have to show that 

F{u,v) < liminf Je,.«(e )(uj ,Wj) . 
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Again, we assume without loss of generality that the sequence ,K(e ) (wj, Wj) 
converges to some finite number c < +00; the claim being trivial if c = +00. 



In particular, this implies that v = 1 almost everywhere. Now (3.8 1, Fatou's 
Lemma, and the F-convergence result for the one-dimensional case imply that, 
for each open set A C 51 and each direction £^ gS^, we have 

liminf J ( A) 

> / F^^y{u^'y,A^^y)dn\y) 



[u- ff + a{£„Vufdx 
+ /?/ H\S,,,.yr\A^^y)d'H\y) 



{u-ff+a{^,Vufdx + p / \{^,Vu)\dV}(x) 



Now let {£,i)ieN C §^ be a dense subset. Then 11, p. 191] implies that 
liminf > ^ / (/ - u)^ + a sup(^j, Vw)^ 



sup|(e„;/„)|rf^""'(y) 



= Fiu) 

3.3 The Upper Bound 



We now prove inequality (3.2 1. To that end, we consider the set yV(r2) consisting 



of all functions u e SBV{Q,) for which the following hold: 

1. Hi(^„\5„) = 0. 

2. The set is the union of a finite number of almost disjoint line segments 
contained in f2. 

3. ^/|o\s„ €lFi^-(fl\5„). 

Obviously, this set is dense in SBV{il) with respect to the L'^-topology. More- 
over, it has been shown in [351 113 that, for every u £ SBV{fl), there exists 
a sequence Uj — >■ u with Uj € W(ri) for every j such that F{uj) — > F{u). 



Using Lemma 3.1 for proving ( |3.2[ ), we therefore have to find for every u e 
>V(il), S > 0, and ej — > sequences uj — > 0, Vj — > 1 as j — >■ 00, such that 
limsupj-^^ >^e,,K(e,)(ui, Wj) < F{u) + 5. 
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Let therefore u € yV(J7) be fixed. By definition of W(il), there exist fc € N 
and ai, hi & 1 < i < k, such that Su = Ui=i[<^i> ^i]- Moreover, 

fc 

n\Su) = n\Su) = J2\\b^-^i\\ ■ 

i=l 

Now define for £ > and < c < 1 

K{e, c) {x e O : dist(a;, Su) < ce} . 

Let moreover /i(e, c) e N and xl^'"-* € Q, 1 < I < fj,{e, c), be such that 

Kie,c)c'[j bM''^^)- 
1=1 

Now note that, if we place the centers of the balls on a line segment [ai,6i], 
then they cover the whole set {.t e : dist(a;, [0^,6^]) < ce} provided that the 
distance between two adjacent centers is at most 2e^/l — . Thus it follows that 
one can cover each set {a; G i7 : dist(a;, [aj, hi\) < ce} with at most ^e^/^z 



balls of radius s. Consequently, we can choose the centers xi in such a way that 
Llie, c) < k + y — — , = fc H , . 



Let now 



1 else. 

Then, for every c we have v^^'"^^ ^ 1 as e — )• 0. Moreover, 
Define now 

' ' [x) := u{x) mm \ ^ . 

Then u'^^''^^ {x) = u{x) for x ^ K{£, c) and u^^^^^^ ?i as e ^ 0. Denoting 
d{x) = dist(a;, S^), we have for almost every x G K{e, c) 

Wu^-'^\x) = Wu{x/-^+u{x)^. 

ce ce 

Thus the triangle inequality and the fact that |Vd(j:)| = 1 almost everywhere 
imply that 

\yu^''-\x)\ < \Vu{x)\^ + < \Vuix)\ + ^ 

ce ce ce 

for almost every x G K{e, c). Therefore, for almost every x G K{e, c), 

\Vu^''^\x)f<2\\Vu\\l + 2^^. (3.10) 
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Now consider each term of the functional Je,K(e) separately. We have 

n Jn\K{e,c) JK{e,c) 



->s^o / {u-ffdx 
Jn 



From (3.10) we get 



/ i;("'^)|Vu("'^)|2dx< / \Vu\^dx+ f K{e)\Vu^'''=^\^dx 

Jn Jn\K(6,c) Jk(s,c) 

Jn\K(e,c) 



Because ^(e) = o(e) as e — > and 

C'^{K{e,c)) < C^{K{e,l)) < 2en\Su) + kne^ = 0{e) as 0, 
it follows that 

limsup / v^^'^^IVw^^'^^prfx < / \\/u\^dx. 
e-i-o Jn Jn\s^ 

Finally, the construction of wf^^^) implies that 



2£m,,„(,)(z;(^^^),r!)<2te + 



Let now ej ^ as j — ^ oo and define Uj :— u^^i^'^\ Vj := v^^^''^\ Then it 
follows that 

limsup Je^.,«(e^)(Mj,Wj) < / {u-ffdx+^ [ |Vup dx+ ^=£=H^(5„) . 
.7-s-oo ^ Jn ^ Jn\s,. Vl - c 



Since < c < 1 was arbitrary and 'H^{Su) < oo, we obtain (3.2) with 6 
(1 — l/Vl — 0^)71^ (Su) , which concludes the proof of Theorem 



3.1 



4 Numerical Implementation 

4.1 Proposed Algorithm 



Based on Theorem [2Tj we propose the following algorithm for the approximate 
minimization of the functional J^ ^ for fixed e > and k > 0. 

Algorithm 4.1 Let f e L°^{i}), a, (3 > 0, e > 0, and < k <1 be given. 
Setk^O and Kg 0. 

1. Define 

\ K if X e Kk, 
Vk(X) := < 

li ifxen\Kk. 
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2. Define 

Uk ■■= argminG(u,Wfc) . 

u 

3. Find y Cz Q\ Kk such that \Vuk{y)\ is maximal. 
I If 

^n\^\Vuk{yr<^, (4.1) 

stop. 

Else set Kk+i '.= Kk U B^{y), increase k by one and go to step 1. 



Steps 3 and 4 of the algorithm use the results of Theorem |2.1[ This theorem 
states that adding in the fc-th step a point y to the edge indicator will approx- 
imately result in a decrease of the functional G by approximately e^a7r(l + 
K)|VMfc(2/)p/(l — k). Thus, we will obtain the steepest descent, if we add a 
point y, where |VMfc(y)| is maximal. At the same time the adding of another 
ball leads to an increase of the term TO£_k by 2/3e. In total, the value of J; 



will increase if (4.1) holds, else « will decrease and therefore it makes sense 



to include the point y into the edge set. 

Remark 4.1 In order to increase the performance of the algorithm, it makes 
sense to add not just one hall in each iteration, hut rather several ones. Also 
in this case a similar approximation as Theorem \2.1\ holds, and thus the same 
criterion for adding new points can he applied. This strategy has been used in 
the numerical examples below. 



4.2 Numerical Results 

We now compare the results obtained with Algorithm |4 . 1 1 with results obtained 
using the approximation introduced by Ambrosio and Tortorelli [3 . This latter 
method consists in minimizing the functional 

I,{u,v) :=l f {u-ffdx+'l f v'\Vu\'dx+^ I [e\Wu\'' + Uv-lf)dx. 

(4.2) 

Again, the function v serves as an edge indicator in the sense that the points 
where v is close to zero are an approximation of the edge set K of the solution 
of the Mumford-Shah functional. In contrast to the approximation by means 
of the functional Je,K, however, where the edge set is given as the points where 



the function v is equal to k, in case of the functional (4.2 1 one has to threshold 
V in order to obtain a precisely defined edge set. 

The minimization of has been carried out by alternately solving the cor- 
responding Euler-Lagrange equations with respect to u and v. For the dis- 
cretization, we have used a finite element approach with piecewise bilinear basis 
functions on the pixel grid. The same discretization has been used for the 



computation of Uk in the second step of Algorithm 4.1 

Figures [T] and [2] show a comparison of the results of the Ambrosio-Tortorelli 
approximation and Algorithm |4.1[ The edge indicators are in both cases com- 



parable, though our algorithm in general classifies more points as edges. The 



22 



main difference between the results is that the Ambrosio-Tortorelh approxima- 
tion leads to a diffuse edge indicator, while our method produces well defined 
edges. As a consequence, also the smoothed images tend to be less blurred; 
compare, for instance, the various light refiections in Figure [2] 
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